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ABSTRACT 

We present 850/im and 450/im data from the JCMT Gould Belt Survey obtained with 
SCUBA-2 and characterise the dust attributes of Class I, Class II and Class III disk sources in 
L 1495. We detect 23% of the sample at both wavelengths, with the detection rate decreasing 
through the Classes from I—III. The median disk mask is 1.6 x 10 -3 M©, and only 7% of Class 
II sources have disk masses larger than 20 Jupiter masses. We detect a higher proportion of 
disks towards sources with stellar hosts of spectral type K than spectral type M. Class II disks 
with single stellar hosts of spectral type K have higher masses than those of spectral type M, 
supporting the hypothesis that higher mass stars have more massive disks. Variations in disk 
masses calculated at the two wavelengths suggests there may be differences in dust opacity 
and/or dust temperature between disks with hosts of spectral types K to those with spectral 
type M. 

Key words: submillimetre: ISM; stars: formation; stars: low mass. 


1 INTRODUCTION 

One of the fundamental properties of circumstellar disks, particu¬ 
larly in relation to planet formation, is their mass (Andrews et al. 
2013). The most useful diagnostic of disk mass is the brightness of 
dust continuum emission at submillimetre wavelengths, where the 
dust emission is optically thin over much of the disk volume, and 
is not overly dependent on temperature. Disk masses are thought 
to decrease during the evolutionary process, either through internal 
processes, or due to interactions with nearby massive stars. 

The traditional classification of pre-main sequence (PMS) 
sources is defined observationally through the Class I—III phases 
(Lada 1987; Adams et al. 1987). Class I sources are more evolved 
embedded young stellar objects (YSOs) with massive accretion 
disks, while Class II sources are T Tauri stars with gas-rich circum¬ 
stellar disks where the supply of envelope material has dissipated. 
Class III sources are stars with tenuous or debris disks, where at 
least the inner part of the disk has been evacuated. The bulk of 
disk mass is accreted during the earliest embedded phases, but more 
evolved sources continue to accrete material onto the disk and stel¬ 
lar host. There is not always a clear distinction between Class I and 
Class II sources, with detailed modelling demonstrating that em¬ 
bedded Class I sources observed face-on can be classified as Class 
II, while Class II disks observed edge-on may look like there is still 
a significant amount of envelope material remaining, and be classi¬ 
fied as Class I (Whitney et al. 2003; Crapsi et al. 2008; Evans et al. 
2009). 


1.1 JCMT GBS 

The Gould Belt Legacy Survey (GBS; Ward-Thompson et al. 2007) 
has been awarded 612 h of time on the James Clerk Maxwell 
Telescope (JCMT) to survey nearby star-forming regions (within 
500 pc), using Heterodyne Array Receiver Programme (HARP; 
Buckle et al. 2009) and Submillimetre Common-User Bolometer 
Array 2 (SCUBA-2; Holland et al. 2013). The GBS will trace the 
very earliest stages of star formation, providing an inventory of all 
the protostellar objects contained in the nearby molecular clouds of 
the Gould Belt, covering ~700 deg 2 . The key science goals of the 
SCUBA-2 survey are: 

• To calculate the duration of each of the protostellar stages. 

• To elucidate the nature of the evolution of protostellar col¬ 
lapse. 


• To discover the origin of the IMF of stars, from intermediate 
mass stars to sub stellar objects. 

• To discern the connection between protostars and the molecu¬ 
lar cloud structure from which they formed. 


1.2 Taurus 

The Taurus star-forming complex is young, with typical YSO ages 
of-1 Myr (Kenyon & Hartmann 1995). As a result, even its most 
evolved sources may have remnants of dusty disks surrounding 
them, since low-mass stars have a median disk lifetime of a few 
Myr (Williams & Cieza 2011). This is one of the key reference 
regions for the distribution of disk masses (Andrews & Williams 
2005; Andrews et al. 2013). It is well-studied, nearby, and contains 
some of the most complete and well-characterised catalogues of 
disk sources (e.g., Furlan et al. 2011; Luhman et al. 2010; Esplin 
et al. 2014). The complex has a significant population of YSOs 
that have no remaining envelope, exposing the disk, and the lack of 
Class 0 sources provides evidence that even the embedded sources 
are reasonably evolved. The observations of L1495 presented here 
probe dust continuum emission which spans protostellar evolution¬ 
ary phases. Moreover, the wide-field mapping shown here ensures 
disk parameters can be examined without selection bias when ex¬ 
amining source statistics, something shown to be problematic for 
previous snapshot surveys (Andrews et al. 2013). 

The low stellar density in Taurus (1-10 stars/pc -2 , Gomez et 
al. 1993; Giidel et al. 2007), with the most massive star in L1495 
an intermediate Ae star, suggests that disk evolution is expected 
to proceed through dynamical star-disk or, in multiple systems, by 
disk-disk interactions, rather than evaporation by radiation from 
nearby stars. Recent observations (Lee et al. 2012; Dodds et al. in 
prep) indicate that the eventual primitive planetary systems that are 
the end results of disk evolution may be dependent on the natal en¬ 
vironment, a topic which can be investigated using the large-scale 
maps obtained by the GBS survey. 

One of the most direct measures of planet-forming capability 
in disks is the mass that they contain. The submillimetre regime, 
where thermal dust emission is optically thin throughout the disk, 
provides one of the best methods of inferring disk mass. At these 
wavelengths, observations are sensitive to micron- to millimetre¬ 
sized grains. With grain growth and dust settling expected to occur 
in the inner disk, where temperatures are warmer, submillimetre 
observations are able to trace the mass in disks to later evolutionary 
phases than infrared observations can (Hernandez et al. 2007; Ricci 
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et al. 2010). With dual-wavelength observations of all the Class I, 
II, and III sources in L1495, we are able to assess the effects of 
temperature and dust opacity on measurements of disk mass in a 
consistent manner. 

Sec. 2 describes the observations, data reduction and process¬ 
ing and disk sample. We present analyses of the physical charac¬ 
teristics of the disks in Sec. 3, and discuss disk evolution in Sec. 4. 
Finally, in Sec. 5 we summarise the conclusions of this research. 


2 OBSERVATIONS 

The SCUBA-2 (Holland et al. 2013) observations presented here 
form part of the JCMT Gould Belt Survey (GBS, Ward-Thompson 
et al. 2007). Continuum observations at 850 /im and 450 jUm 
were made using fully sampled 30' diameter circular regions 
(PONG1800 mapping mode; Chapin et al. 2013). Also included 
are observations from the GBS Science verification phase, taken in 
PONG900 and PONG3600 mapping modes. Larger regions were 
mosaicked with overlapping scans. 

The data were reduced using an iterative map-making tech¬ 
nique (makemap in SMURF; Chapin et al. 2013), and gridded to 
6" pixels at 850 jUm and 4" pixels at 450 gm. The iterations were 
halted when the map pixels, on average, changed by <0.1% of the 
estimated map rms. The initial reductions of each individual scan 
were coadded to form a mosaic from which a signal-to-noise mask 
was produced for each region. The final mosaic was produced from 
a second reduction using this mask to define areas of emission, with 
individual scans cropped to exclude noisy edge regions. 

Detection of emission structure smaller in size than the spatial 
filter are robust both inside and outside of the masked regions. A 
spatial filter of 600" is used in the reduction, and tests with simu¬ 
lated data show that flux recovery is robust for sources with a Gaus¬ 
sian FWHM less than 2.5', provided that the source peak brightness 
is several times the noise. Sources between 2.5' and 7.5' will be de¬ 
tected, but both the flux and the size are underestimated because 
Fourier components with scales greater than 5' are removed by the 
filtering process. Detection of sources larger than 7.5' is dependent 
on the mask used for reduction. In this analysis, we focus on com¬ 
pact disk sources, which, even with their remnant envelopes, have 
size scales <1' (0.04 pc at a distance of 140 pc), where we expect 
flux recovery to be robust at both 850 gm and 450 gm. 

The data are calibrated in mJy pixel -1 , using aperture FCFs 
of 2.34 Jy/pW/arcsec 2 and 4.71 Jy/pW/arcsec 2 at 850 gm and 
450 gm, respectively, derived from average values of JCMT cal¬ 
ibrators (Dempsey et al. 2013), and correcting for the pixel area. 
The PONG scan pattern leads to lower noise in the map centre and 
overlap regions, while data reduction and emission artefacts can 
lead to small variations in the noise over the whole map. 

Fig. 1 shows the calibrated emission maps at 850 gm and 
450 gm on the left, and the error maps on the right, which are also 
contoured with the mask used for data reduction. A mean 1 a un¬ 
certainty of 1.7 mJy pixel -1 at 850 gim was achieved in the final 
L 1495 mosaic. The inclusion of science verification data leads to 
a variation in noise across the final mosaic, due to the differences 
in spatial coverage. We achieve a better, more consistent sensitiv¬ 
ity within the masked regions, however. Within the masked region, 
which is the focus of the science presented here, a mean 1 a rms of 
0.98 mJy pixel -1 , equivalent to 2.5 mJy beam -1 , was achieved. For 
the 450 gm data, a mean 1 a rms of 5.1 mJy pixel -1 , equivalent to 
13.3 mJy beam -1 , was achieved inside the masked region. 


2.1 Disk sources in L1495 

We use the most comprehensive disk catalogues of Taurus (Luhman 
et al. 2010; Esplin et al. 2014) to identify Class I, Class II, and Class 
III sources in L1495, which are listed in Table 1. We have applied 
updated spectral types published by Herczeg & Hillenbrand (2014) 
and Pecaut & Mamajek (2013), as well as details of multiplicity of 
PMS stars in Taurus (e.g., Akeson & Jensen 2014) and the classifi¬ 
cation of sources as transition disks (e.g., Najita et al. 2007; Furlan 
et al. 2011). Note that there are no known Class 0 sources in this re¬ 
gion. Figures A1-A56 in Appendix A show thumbnail continuum 
images of all the disks in our sample. 


2.2 CO contamination of dust continuum 

Dust continuum measurements in the submillimetre, observed us¬ 
ing wide-band filters, can be contaminated by flux from molecular 
lines. Lines of 12 CO are particularly problematic, as they are the 
brightest lines in the submillimetre wavelength range. For SCUBA- 
2, the 12 CO J = 3 2 transition is near the central transmis¬ 

sion peak of the 850 gm bandpass filter, and can dominate dust 
continuum observations in regions of energetic molecular outflow 
(Drabek et al. 2012). The JCMT GBS datasets include large-scale 
molecular line maps in 12 CO J = 3—^2, which we utilise, follow¬ 
ing the techniques described in Drabek et al. (2012) to quantify 
the contribution of 12 CO emission to our dust continuum measure¬ 
ments. 

We use the JCMT GBS HARP observations of Taurus (Davis 
et al. 2010), which comprise CO J = 3 —)► 2 maps along the north¬ 
west region of Taurus known as L1495. The region contains a dif¬ 
fuse cloud (named as a bowl or hub, Davis et al. 2010) which is 
the focus of the dust continuum observations presented here. The 
CO data were processed through the same iterative mapmaker as 
the 850 gm data, to apply the same spatial filtering (Hatchell et al. 
2013; Sadavoy et al. 2013). 

In the quiescent Taurus region, emission from CO is relatively 
diffuse, with only a few locations of strong emission, associated 
with high velocity protostellar outflows (Davis et al. 2010). The 
level of CO contamination does not exceed 15%, the 850 gm sys¬ 
tematic flux uncertainty, and is generally below 5%, as shown in 
Table 1 for the disk sources in our sample. Sources not detected at 
850 gm are marked as while sources which were not observed 
in 12 CO are marked as ‘N/0’ in the Table. The results are consistent 
with results from other clouds (e.g. Drabek et al. 2012), where high 
levels of CO contamination in the 850 gm dust continuum emission 
are associated with energetic molecular outflows. Such outflows are 
rarer in the L1495 region observed here, which is more evolved 
than other parts of Taurus (Davis et al. 2010). Of the two Class I 
and three Class II disk sources with CO contamination levels >5%, 
IRAS 04111+2800 G (Class I) drives a high velocity outflow, while 
CoKu Tau 1 (Class I), MH03 (Class II) and V410 X-ray 6 (Class II) 
are associated with high velocity red-shifted CO emission (Davis et 
al. 2010). The remaining Class II source with contamination >5%, 
LR1, is a compact dust source embedded in weaker extended emis¬ 
sion (both dust emission and CO integrated intensity) with local 
CO spectra having multiple velocity components. 

Given the low levels of CO contamination in the 850 gm band, 
we do not expect any measurable contamination in the 450 gm 
band, where the 12 CO / = 6 —)► 5 transition lies at the edge of the 
SCUBA-2 bandpass filter (Drabek et al. 2012). 
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Right ascension (J2000) Right ascension (J2000) 


Figure 1 . SCUBA-2 850 ]Um (top) and 450 /im (bottom) continuum images (left) and error maps (right) of L1495. Fluxes are in mJy pixel -1 . Contours on the 
error maps show the mask used in the data reduction (Sec. 2). Note that the data maps are on a logarithmic colour scale, while the error maps are on a linear 
scale, with the colour-coding reversed . 


3 ANALYSIS 

Background subtracted mean and integrated flux densities for the 
disk sample were measured using 20" diameter apertures, with the 
background level estimated from nearby pixels in annuli between 
30" and 40". For a few crowded regions, nearby background re¬ 
gions with 20" diameter apertures were chosen manually. Table 1 
shows the integrated flux values with 1 a uncertainties for each disk 
source, as well as information on the Class and spectral type. 

Using a detection criteria of 3 <7 or better in total flux, we 
detect 88% (7/8) of the Class I sources at each wavelength. For 
Class II sources, we detect 38% (11/29) at 850 /n m and 21% (6/29) 
at 450 jUm. Only 11% (2/19) of the Class III sources are detected 
850 gm, and none at 450 jUm. Table 2 summarises detections at the 
3 a level for all the sources, by wavelength, Class, spectral type and 
multiplicity. As discussed in Sec. 2, our reduction algorithm filters 
out emission on scales of a few arcminutes, and so these observa¬ 
tions are not sensitive to emission from the surrounding cloud. The 
presence of large remnant envelopes can increase the submillimetre 
continuum around disk sources, particularly for more embedded, 
younger sources. Removing the local background using close an¬ 


nuli should compensate in part for such contributions, although the 
results are then more likely to be affected by contamination from 
the error beam. We investigate the contribution of emission from 
surrounding envelopes and loud material in Sec. 3.2. 

The sources in our sample are dominated by disks with stellar 
hosts of spectral type M, which have lower detection rates than for 
disks with hosts of spectral type K, as previously found in stud¬ 
ies focussing on late spectral types (Andrews et al. 2013; Schaefer 
et al. 2009). Several of the sources are known binary or multiple 
sources, and these have a slightly higher detection rates (29% and 
22% for multiple and single sources, respectively, at both wave¬ 
lengths) than sources thought to be single systems, although there 
is, again, a large difference in the numbers of sources in each cat¬ 
egory. There are no known multiple Class I sources in our sample, 
and we do not detect the single Class III multiple source. We anal¬ 
yse these trends further in Sec. 4. For sources that are not detected, 
a flux limit can be obtained from the mean and standard error on 
the mean of the associated integrated flux values, listed in Table 
1. For Class II and III sources at 850 jUm, the mean flux and stan¬ 
dard error on the mean are 6.24 ± 2.38 mJy and 1.55 db 1.57 mJy, 
respectively. At 450 jUm, the mean flux and standard error on the 
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Table 1. Background-subtracted continuum flux measurements, from aperture photometry, with 1 a uncertainties on the mean and integrated fluxes, which do 
not include the calibration errors. Sources are ordered by Class, then Right Ascension, with names, Class and spectral types as given in the catalogues (Luhman 
et al. 2010; Esplin et al. 2014) 


#a 

ID 

Class 

Sp.Type 

mean F 850 
mJy pixel -1 

total F 85 o 
mJy 

mean F 450 
mJy pixel -1 

total F 450 
mJy 

%CO 

1 

J04153566+2847417 

I 

- 

2.85±1.17 

25.7±10.6 

4.25±7.16 

82.9±140 

N/O 

2 

SSTTau04183 lp2+28261 

I 

- 

3.40±0.36 

30.7±3.27 

6.20±1.32 

120±26 

<5 

3 

IRAS04108+2803B 

I 

K0 

5.91±0.44 

53.3±4.01 

14.8±3.14 

288±61.6 

<5 

4 

IRAS04111+2800G 

I 

K0 

4.60±0.50 

41.5±4.54 

11.0±2.94 

215±57.8 

13 

5 

IRAS04108+2803A 

I 

M4.75 

4.63±0.44 

41.8±3.99 

14.7±3.14 

287±61.6 

<5 

6 

IRAS04154+2823 

I 

M2.5 

7.20±0.43 

65.0±3.89 

ll.3dzl.37 

221±26.9 

N/O 

7 

CoKuTau/1 

I 

M0 

1.91±0.40 

17.3±3.64 

4.58±1.25 

89.3±24.6 

8 

8 

IRAS04158p2805 

I 

M5.25 

22.3±0.46 

201±4.15 

56.7±2.24 

1110±44.0 

<5 

9 

V410X-ray7 

II 

M0.75 

1.38±0.36 

12.5±3.27 

2.61d=1.37 

50.9±27.0 

<5 

10 

J04141188+2811535 

II 

M6.25 

0.94±0.44 

8.50±3.96 

-3.84±3.13 

-74.9±61.6 

<5 

11 

V773 Tau A+B 

II T 

K3 

1.23±0.44 

11.1±3.96 

-6.02±3.14 

-117±61.7 

<5 

12 

FMTau 

II 

M4.5 

1.19±0.48 

10.8±4.37 

-2.17±2.88 

-42.3±56.5 

<5 

13 

FNTau 

II 

M3.5 

1.62±0.53 

14.7±4.78 

8.07±3.11 

158±61.0 

N/O 

14 

CWTau 

II 

K3.0 

6.55±0.49 

59.2±4.41 

16.9±2.90 

330±57.0 

<5 

15 

CIDA1 

II 

M4.5 

0.39±0.50 

3.56±4.58 

4.59d=2.99 

89.6±58.6 

<5 

16 

MHO 1/2 

II 

M2.5 

48.9±0.46 

442±4.20 

85.7±3.03 

1670±59.5 

<5 

17 

MH03 

II 

M2.2 

2.02±0.51 

18.3±4.62 

2.30±3.03 

44.8±59.6 

10 

18 

FO Tau A+B 

II T 

M3.9 

1.80±0.46 

16.2±4.20 

7.82±2.73 

153±53.5 

N/O 

19 

J04153916+2818586 

II 

M4.0 

1.64±0.78 

14.8±7.04 

1.36±5.59 

26.5±110 

<5 

20 

J04155799+2746175 

II 

M5.2 

2.31±2.02 

20.8±18.3 

13.9±11.1 

271±217 

N/O 

21 

J04161210+2756385 

II 

M4.75 

-0.28±0.84 

-2.54±7.66 

-6.06±6.04 

-118±119 

N/O 

22 

J04163911+2858491 

II 

M5.5 

-0.33±1.90 

-3.02±17.20 

3.82±7.01 

74.6±137 

N/O 

23 

CYTau 

II 

M2.3 

13.7±0.30 

124±2.68 

10.0dil.38 

196±27.0 

<5 

24 

KPNO10 

II 

M5 

2.21±0.44 

20.0±4.03 

-1.35±1.94 

-26.3±38.1 

<5 

25 

Y41 OX-ray 1 

II 

M3.7 

-0.12±0.45 

-1.06±4.09 

-0.03dzl.99 

-0.60±39.0 

N/O 

26 

V410Anonl3 

II 

M5.75 

-0.68±0.41 

-6.16±3.75 

-0.01±1.31 

-0.28±25.7 

N/O 

27 

DD Tau A+B 

II 

M4.8 

2.10±0.38 

19.0±3.41 

6.17±1.43 

121d=28.1 

<5 

28 

CZ Tau A+B 

II 

M4.2 

-0.06±0.37 

-0.56±3.38 

-0.77±1.43 

-15.0d=28.0 

<5 

29 

Y410 X-Ray 2 

II 

M0 

1.25±0.42 

11.3±3.84 

0.59d=1.31 

11.6±25.8 

N/O 

30 

V892Tau 

II 

A0 

43.3±0.34 

391±3.06 

74.3±1.24 

1450±24.3 

<5 

31 

LR 1 

II 

K4.5 

3.05±0.40 

27.5±3.65 

ll.7dzl.23 

229±24.1 

13 

32 

V410X-Ray 6 

II T 

M5.9 

1.10±0.43 

9.90±3.94 

3.53±1.44 

68.9±28.3 

9 

33 

KPN012 

II 

M9 

0.33±0.97 

2.96±8.77 

4.25±6.18 

82.9±121 

N/O 

34 

FQ Tau A+B 

II T 

M4.3 

-0.31d=0.74 

-2.77±6.73 

0.70±5.60 

13.6±110 

N/O 

35 

V819Tau 

II T 

K8.0 

-0.33±0.75 

-2.94±6.83 

10.4±5.68 

203d=112 

N/O 

36 

FRTau 

II 

M5.3 

0.25±0.77 

2.26±6.98 

4.22±5.59 

82.4±110 

N/O 

37 

J04201611+2821325 

II 

M6.5 

3.83d=2.33 

34.6±21.1 

2.39d=12.1 

46.6±238 

N/O 

38 

LkCal 

III 

M4 

-0.26±0.51 

-2.38±4.64 

-2.70±3.10 

-52.7±60.8 

<5 

39 

Anonl 

III 

M0 

-0.64d=0.48 

-5.74±4.37 

0.90±2.93 

17.5±57.5 

<5 

40 

XEST20-066 

III 

M5.2 

0.42±0.54 

3.77±4.87 

-4.96±3.01 

-96.8±59.0 

<5 

41 

XEST20-071 

III 

M3.1 

1.34±0.48 

12.1±4.38 

0.02±2.81 

0.42±55.1 

<5 

42 

CIDA2 

III 

M5.5 

0.63±0.48 

5.71±4.34 

-0.33±2.79 

-6.35±54.7 

<5 

43 

KPNOl 

III 

M8.5 

0.10±0.91 

0.93±8.27 

0.53±5.37 

10.3d=106 

N/O 

44 

J04161885+2752155 

III 

M6.25 

-0.79±0.82 

-7.17±7.43 

3.89±5.93 

75.9±116.3 

N/O 

45 

LkCa4 

III 

Ml.3 

1.03±0.81 

9.32d:7.37 

16.0±6.21 

312±122 

N/O 

46 

LkCa5 

III 

M2.2 

-0.54±0.48 

-4.91±4.32 

-3.85±2.10 

-75.2±41.3 

N/O 

47 

V410Xray3 

III 

M6.5 

-0.87±0.41 

-7.88±3.69 

-1.74±1.30 

-33.9±25.6 

<5 

48 

V410Anon25 

III 

Ml 

1.68±0.36 

15.18±3.26 

1.81±1.32 

35.4±25.9 

<5 

49 

KPNOll 

III 

M5.9 

-0.08±0.47 

-0.68±4.28 

0.54±2.23 

10.5±43.8 

N/O 

50 

V410Tau A+B+C 

III 

K7 

0.43±0.40 

3.91±3.59 

-0.24±1.24 

-4.76±24.3 

<5 

51 

V410X-ray4 

III 

M4 

1.94±0.39 

17.50±3.58 

-0.98±1.21 

-19.2±23.8 

<5 

52 

Hubble4 

III 

K8.5 

-0.17±0.40 

-1.52d=3.59 

0.61±1.24 

12.0±24.3 

<5 

53 

KPN02 

III 

M7.5 

0.49±0.46 

4.46d=4.21 

-4.92±1.94 

-96.0±38.1 

<5 

54 

V410X-ray5a 

III 

M5.5 

-0.17±0.43 

-1.50±3.95 

-0.93dzl.41 

-18.1±27.7 

N/O 

55 

J04144158+2809583 

III 

- 

0.94±0.48 

8.45±4.33 

0.56±2.81 

10.9±55.1 

<5 

56 

J04161726+2817128 

III 

- 

1.26±0.80 

ll.4dz7.28 

0.78±6.03 

15.1±118 

<5 


a 850 /im detected sources are indicated in boldface, and underlined if also detected at 450 /im 
T identifies a transition disk 
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Table 2. 3 a detection statistics in percentages for the disk sources. 


Category # sources % Detections 

850 gm 450 pm 


All 

56 

36 

23 

I 

8 

88 

88 

II 

29 

38 

21 

III 

19 

11 

0 

Sp Type K 

8 

50 

50 

Sp Type M 

43 

33 

16 

Multiple 

7 

43 

29 

Single 

49 

35 

22 


mean are 42.7 d= 19.7 mJy and 5.1 ± 19.8 mJy for Class II and III 
respectively. 


3.1 Physical characteristics of dust in the disk sample 

A measure of the planet-forming capabilities of circumstellar disks 
is determining if they contain sufficient mass (>20Mjup of dust 
and gas) to create the equivalent of a primitive Solar System (Davis 
2005). Accurately measuring disk masses is a vital step in under¬ 
standing disk evolution and the diversity of exoplanets. The opacity 
of dust is known to be one of the largest uncertainties in accurately 
assessing disk mass (e.g. Beckwith et al. 1990). With observations 
at two wavelengths for our complete sample of disks we investigate 
the spectral slope, and through it the dust opacity and temperature, 
before calculating disk masses. The analyses are described in the 
following sections, and the results presented in Table 3. Results and 
calculations are all for detections at >3 <7 unless otherwise stated. 


3.1.1 Continuum slopes 


The slope ( a ) of the continuum submillimetre emission can be cal¬ 
culated from the simple spectral index: 

ln(F 450 ) -ln(F 850 ) 

ln(850) — ln(450) V ’ 

for sources with emission detected at 3 cr or greater in both wave¬ 
lengths. For the disk sources, we use the integrated fluxes re¬ 
ported in Table 1 to calculate individual a values. For optically thin 
and isothermal emission, where the Rayleigh-Jeans criteria applies 
(e.g., dust temperatures >20 K for hv >> kT at 850 /im, and >35 K 
at 450 /im), then a=2+p, where j8 is the power law index of the 
dust opacity (k) at submillimetre wavelengths: 


k*v i — k v 2 



( 2 ) 


In this case, we expect a values of 3 if jS=1.0, as appropriate for 
disk sources (e.g. Beckwith et al. 1990), or tending to higher val¬ 
ues if P >1.0. For optically thick emission, or for dust grains much 
larger than the wavelength, values of a ^2 are expected, regardless 
of P values. To ensure the same scales are sampled at both wave¬ 
lengths, we use cross-convolved datasets for this, by convolving the 
850 gm map with 2-dimensional Gaussians representing the pri¬ 
mary and secondary beams at 450 gm, and similarly for the 450 pm 
map to produce resolution-matched maps. The FWHM of the pri¬ 
mary/secondary beams are 7.9 // /13.0 // and 25.0 // /48.0 // at 450 pm 
and 850 /im, respectively. The error beam at 450 pm contributes 



Figure 2. Normalised distribution of a values (left axis, stepped histogram) 
calculated from the simple spectral index (Eq. 1) for all pixels detected with 
>3 a in the L 1495 region. Individual counts (right axis) for Class I (grey) 
and Class II (white) disk sources are also shown. 


significantly to flux in a large area (Dempsey et al. 2013), where the 
relative volumes of the primary and secondary beams are 0.78/0.22 
at 850 pm and 0.61/0.39 at 450 gm. Fig. 2, (left axis, stepped his¬ 
togram), shows the distribution of a across L1495, obtained from 
all pixels in Fig. 1 with a detection >3 a. The distribution has a 
mean and standard deviation of 2.8±0.6. Since regions with sig¬ 
nificant detections at 450 pm are compact, and generally towards 
star-forming cores in the L 1495 region, this distribution suggests 
we are tracing little of the extended cloud material, where P values 
of 1.5 (and a ^3.5 for a=2+P) or higher may be expected. 

Fig. 2 also shows the distributions of a values in the disk 
sources, which are similar for Class I and Class II sources. CY Tau 
has an anomalously low value of 0.42, although this is consistent 
with previously published submillimetre values (e.g., Andrews & 
Williams 2005). There are four sources with a > 2.5, implying 
P ~1.0, and eight sources with a < 2.5, implying either grain 
growth, or optically thick emission. Although we are not able to 
distinguish between these scenarios, the similar distribution we find 
between Class I and Class II disks provides evidence that disk 
evolution must be rapid. Under the assumption of optically thin, 
isothermal and warm (>20K) emission from dust in the disks, 
these results suggest that a value of P=l.O is appropriate to use 
in the rest of the analyses presented here. 


3.1.2 Temperature 

Alternatively, the dust temperature can be inferred from the flux 
ratio (Eq. 3), given an assumed value of /3: 

^50 = /850 \ 3+/3 / exp(hcA 85 ofer dust )-l \ 

^850 V 450 J V exp (he/ A 45 o&7d ust ) - 1 ) 

A flux ratio map was created using the cross-convolved datasets, as 
described in Sec. 3.1.1 The dust temperature map was created from 
the ratio map for /3=1.0, as determined in Sec. 3.1.1, with a high 
temperature cap of 50 K, above which the flux ratio does not con¬ 
strain the dust temperature (e.g., Salji et al. 2014). Fig. 3 (left axis) 
shows the normalised temperature distribution calculated from the 
flux ratio (Eq. 3) of the L1495 region for P=l.O and /3=1.8. A value 
of /3=1.8 is more appropriate for emission from the ISM, while 
the disk sources are more accurately described with P=1.0. Also 
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Figure 3. Normalised distribution of dust temperatures (left axis) calculated 
from the flux ratio map (Eq. 3), assuming /3=1.8 (blue) and j8=1.0 (black). 
Temperature counts for Class I (grey) and Class II (white) disk sources are 
shown on the right axis, assuming /3=1.0. 



Figure 4. Mass distribution of disk sources, assuming T dust =20 K. Class I, 
II and III are shown in grey, white and black respectively. One Jupiter mass, 
Mj UP = 1.05 x 1O -3 M 0 . 


shown is the temperature distribution for Class I and Class II disk 
sources, calculated assuming /3=1.0. Although both Classes have 
similar distributions of temperatures, the sources with the high¬ 
est temperatures are all Class I. The median flux ratio for all the 
sources is 4.3, implying a dust temperature of 20 K if jS=1.0. The 
Class I sources’ flux ratios extend to higher values than those of 
the Class II sources, with median flux ratios of 4.6/4.0 leading to 
dust temperatures of 23/18 K, respectively. We are removing the lo¬ 
cal background emission in which the disk sources are embedded. 
However, it is possible that the envelopes may be still contributing 
to flux from the Class I sources. If so, this will affect the flux ratios 
to a larger extent than for the Class II sources. Higher 450/850 flux 
ratios can be due to higher dust temperatures and/or a higher opac¬ 
ity index j8 (Hatchell et al. 2013), and requires flux measurements 
at additional wavelengths to disentangle. 

Two disk sources have 5 K^T dust ^10K, while three have 
30K^T ( j ust <50K, and one has a flux ratio too high to estimate 
the temperature reliably. Note, however, that the majority of the 
disk sources are in the regime where the Rayleigh-Jeans criteria 
fails at one or both wavelengths, and so the continuum slope is not 
well constrained. These sources may have contributions from sur¬ 
rounding material with different dust emission properties, requiring 
different /3 values, and for these sources interferometric observa¬ 
tions will be required to isolate emission from the compact disk. 
The similar range of dust temperatures calculated for both Classes 
suggests that adopting a single constant dust temperature for all 
sources not detected at both wavelengths will provide consistent 
results. A dust temperature of 20 K is consistent with the median 
values derived for both Class I and Class II sources, and also in 
line with previous work (e.g., Mohanty et al. 2013; Andrews & 
Williams 2005; Williams et al. 2013). We therefore adopt a con¬ 
stant disk temperature of 20 K in calculating disk masses for the 
rest of the analyses presented here. 


3.1.3 Masses 

To determine the disk mass for all of the sources detected in at least 
one wavelength, we use values of ^^=20 K and /3=1.0, consistent 
with the above analyses. We assume that the emission is optically 
thin, the opacity law does not change from source to source, and 


can be described with a single characteristic temperature across the 
unresolved disk, so that the disk mass (dust + gas) scales with the 
observed flux density (. F v ): 


M d = 


F V D 2 

K\rB\< (7j U st) 


(4) 


where D is the distance to Taurus (we adopt the canonical value 
of 140pc, Elias 1978). B v (T dust ) is the Planck function for the 
given dust temperature, and K v is the wavelength-dependent dust 
opacity. We adopt a dust opacity of ^850 = 0.0035 m 2 kg -1 , which 
assumes a gas-to-dust ratio of 100:1, and follows from the com¬ 
monly adopted value for disks of Ki.3mm = 0.0023 m 2 kg 1 (e.g., 
Beckwith et al. 1990), with /3 = 1 (Eq. 2). Where sources were de¬ 
tected at both wavelengths, masses have also been calculated from 
the inferred dust temperature (Eq. 3), and these are both listed in 
Table 3. 

Fig. 4 shows the mass distribution of the disk sources, assum¬ 
ing T dust =20K, which has a mean disk mass of 4.5 xl0 _3 M Q 
(~4Mjup), and a median disk mass of 1.6 x 10 -3 M 0 (~1 Mjup). 
The Class II sources have higher mass, on average, than the Class 
I sources, dominated by the bright A0 source, V892Tau, and the 
multiple source MHO 1/2, both with disk masses > 20 Mjup- 

For sources where we are able to derive a dust colour 
temperature, the two cold sources, with T dust <10K, have 
M(T dust ) > 20 Mjup (assuming j6 = 1). One of these, CYTau, has 
been shown in recent high resolution molecular line observations to 
have a compact disk, at a low inclination (Guilloteau et al. 2014). 
If the disk is compact, then the large mass in a small disk indicates 
high column densities, and so the assumption of optically thin dust 
emission may not be valid, and the disk mass underestimated. 

The above analysis demonstrates that the simultaneous 
dual wavelength wide-area mapping observations attainable with 
SCUBA-2 provide valuable constraints on the physical characteris¬ 
tics of the dusty circumstellar disks. To further demonstrate the va¬ 
lidity of these results, we can compare disk masses derived above 
with those derived using SED models of young stellar objects. 
We use the Web-based tools of Robitaille et al. (2007), obtain¬ 
ing shorter wavelength photometric data from 2MASS, IRAC and 
MIPS using published catalogues (Cutri et al. 2013; Redbull et al. 
2010; Roeser et al. 2008; Cutri et al. 2003). The SED-fitting tool 
(Robitaille et al. 2007) utilises a large grid of pre-computed mod- 
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Figure 5. The filled circles show the input fluxes and the coloured lines 
show the best fits to the SED (Robitaille et al. 2007). The dashed line shows 
the stellar photosphere corresponding to the central source of the best fitting 
model, as it would look in the absence of circumstellar dust (but including 
interstellar extinction). Disk masses derived from the best-fit models of the 
SED are 7.34x 1O - 3 M 0 (CWTau) and 27.9 xlO" 3 M 0 (CYTau). 


Table 3. Physical parameters derived from continuum measurements 


# ID 

Class 

Mass (20K) 
x 10 _3 M o 

Tdust 

K 

Mass (Td us t) 
x 10 _3 M o 

a 

2 

I 

1.69 

24.5 

1.27 

2.69 

3 

I 

2.94 

48.9 

0.91 

2.36 

4 

I 

2.29 

37.6 

0.98 

1.80 

5 

I 

2.31 

50.0 

0.70 

2.91 

6 

I 

3.59 

6.96 

28.0 

1.38 

7 

I 

0.95 

16.8 

1.24 

2.92 

8 

I 

11.1 

14.6 

18.3 

2.32 

9 

II 

0.69 

5.47 

- 

- 

13 

II 

1.05 

16.1 

- 

- 

14 

II 

3.26 

18.2 

3.78 

2.47 

16 

II 

24.4 

17.3 

30.3 

1.73 

17 

II 

1.01 

- 

- 

- 

18 

II 

0.86 

- 

- 

- 

23 

II 

6.84 

7.55 

43.3 

0.43 

24 

II 

1.10 

- 

- 

- 

27 

II 

1.05 

16.7 

1.38 

2.32 

30 

II 

21.6 

12.8 

44.9 

1.80 

31 

II 

1.52 

32.2 

0.79 

2.59 

48 

III 

0.84 

- 

- 

- 

51 

III 

0.97 

- 

- 

- 


els that span a wide range of parameter space in order to fit the 
SED, and provides the resulting model fits as datafiles and as SED 
plots, shown in Fig. 5 for two of the sources detected at both wave¬ 
lengths, CY Tau and CW Tau. The SCUBA-2 observations suggest 
that CY Tau has a more massive disk than CW Tau, which results 
from SED fitting support. The masses derived from SCUBA-2 ob¬ 
servations, and those derived from SED fitting, agree to within a 
factor of 2 for both sources. 

We are also able to calculate masses from the 450 gm flux, 
which should be similar to the masses derived from 850 gm fluxes, 
if our assumptions on the dust opacity index (Eq. 2) are correct 
and extend to 450 gm. For a constant dust temperature of 20 K, the 
mean and standard error of the M 450 /M 850 ratio is 1.16±0.11. The 
lower mass disks have the largest variation in the mass ratio, but 
differences in derived dust temperature with respect to the char¬ 
acteristic temperature of 20 K do not account for this difference. 
The mass ratios with the largest deviations from 1.0 are for those 
sources with a<2 or a>3, and indicate that the largest uncertain¬ 
ties in the mass estimates are from variations in the dust opacity. 
We investigate variations in the mass ratio with respect to source 
and stellar host characteristics in Sec. 4. 


Almost two-thirds (63%) of the sources were not detected 
at the 3 a level or above. The mean and standard error of mass 
for these sources, calculated using the mean of fluxes for sources 
not detected at the 3 a level, reported in Table 1, and using a dust 
temperature of 20K, is 1.07±0.12 x 10 -3 M 0 . We have a similar 
number of non-detections for the Class II and Class III sources, 
which have mean and standard errors of 1.26±0.21 xlO _ 3 M 0 
and 0.83±0.06 xlO _ 3 M 0 . We investigate further the mass lim¬ 
its and structure of non-detections in Sec. 3.2. Of the five Class II 
sources characterised as probable transition disks, we detect one, 
the multiple source FOTau A & B, and calculate a disk mass 
of 0.86 xlO 3 M 0 . Multiplicity in disk sources hampers efforts 
to characterise their physical characteristics (Najita et al. 2007), 
which we investigate in Sec. 4. 

3.2 Source characteristics by Class 

The detection rate for our disk sample decreases with evolution¬ 
ary Class. For example, at 850 gm, there is only one non-detection 
for the eight Class I sources, but only two detections for the nine¬ 
teen Class III sources. The wide-field mapping available with the 
SCUBA-2 datasets means that we can investigate in a statistical 
manner the material surrounding each of the disk sources. By stack¬ 
ing the data at the nominal disk positions, we can compare average 
properties by Class. 

Figs. 6-8 show the stacked images at 850 gm (left) and 
450 gm (right), with no local background subtracted. Also shown 
are the radially averaged emission profiles over the central 2' of 
the stacked image and the relative contributions of the primary and 
secondary (error) beams at each wavelength. 

Fig. 6 shows the stacked images and radially-averaged emis¬ 
sion profiles for the detected Class I sources. At both wavelengths, 
we see compact emission peaks within extended low-level emis¬ 
sion. The apertures used to measure background-subtracted fluxes 
for these disk sources, with a diameter of 20 " and local background 
measured in annuli of 30 // -40 // , probe size scales of 2800 AU- 
5600 AU, and should reduce contamination in the observed disk 
flux due to emission from the envelope, which could extend to 
scales —20000 AU. A 2-component Gaussian fit to the radial pro¬ 
files indicates that the low-level wide component - associated with 
the envelope - will contribute 30% to the fluxes at both wavelengths 
if not removed. The compact peak can be associated with an unre¬ 
solved disk, with widths of 6.3"-7.7" (-1000 AU). At SCUBA-2’s 
resolutions, we are not able to resolve any accretion disc that might 
be present. 

For the Class II sources, we can compare the structure of de¬ 
tected and non-detected sources using stacked data, shown in Fig. 7. 
Nearby compact sources are bright enough to be seen in the Class 
II stacked data for detections and non-detections, but are clearly 
separated from the nominal disk positions. For the detected and 
non-detected sources at both wavelengths, the compact peak has a 
width of 6.3"-7.0" (—1000 AU). At 450 gm, the stacked data for 
non-detections show a weak central peak, which is offset by just a 
few arcseconds from the nominal central position, making it diffi¬ 
cult to determine the contribution to the measured disk flux. This 
offset may indicate dust asymmetry in the disk, suggesting that the 
non-detected Class II disks are more evolved, rather than just less 
luminous. At 850 gm, and for detected 450 gm sources, the narrow 
component contributes at least 95% of the flux. 

The stacked data for the Class III sources are dominated by 
extended arcs of emission at offsets of — 1 ' from the nominal cen¬ 
tral position, which arise in a small number of sources located near 
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Figure 6 . Top: Stacked image of detected Class I sources at 850 pm (left) 
and 450 pm (right), with flux shown in mJy pixel -1 . Bottom: Radially av¬ 
eraged emission profiles for the stacked images, showing the low-level ex¬ 
tended emission in which the compact disk is embedded. Also shown are 
the relative contributions of the primary and secondary (error) beams at each 
wavelength. Using a 2-component Gaussian fit to model the radial profiles, 
the extended emission will contribute 30% to the observed flux if not re¬ 
moved through background subtracted aperture photometry. 



the strongest emission in the map, the bright reflection nebula illu¬ 
minated by the Class II Ae star V892 Tau. These sources (V410 Tau 
and V410 X-Ray 3 undetected at both wavelengths, V410 X-Ray 4 
and Y410 Anon 25 undetected at 450 pm) have been removed from 
the final stacked data, shown in Fig. 8. Extended emission at dis¬ 
tances >1' from the disk position remains in the stacked data, how¬ 
ever. The emission decreases radially towards the disk position, and 
at offsets <0.5' from the disk position, the flux is more consistent 
with a flat profile, particularly at 850 pm. There may be emission 
from the central 6"-10" region, but higher sensitivity data are re¬ 
quired to investigate whether or not the dust ‘hole’ for these sources 
is real. 

With the disks in both Class I and Class II sources unresolved 
at this resolution, we are unable to determine differences in observ¬ 
able disk size at submillimetre wavelengths between the Classes. 
The differences in emission structure that we do observe are instead 
dominated by envelope and/or cloud material. The derived physi¬ 
cal characteristics show only small differences on average between 
Class I and Class II objects, with those in each phase covering sim¬ 
ilar ranges in temperature, a , and disk mass. In the next Section, 
we investigate alternative indicators of disk evolution. 


4 DISK EVOLUTION 

Disk masses are known to decline with age, and to be intrinsically 
lower for late-type stellar hosts. Andrews et al. (2013) have shown 
that the disk luminosity and mass are correlated with stellar host 
mass for Class II sources, with more massive stars more likely to 
have more massive disks. In this section, we investigate this corre¬ 
lation for M- and K- spectral type Class I and Class II sources in 
L1495. 

A wide variety of stellar masses have been computed from 
spectral types of young, pre-main sequence stars, and to date it is 
been difficult to determine dynamical stellar masses (Schaefer et 
al. 2009). We therefore use stellar types as a proxy for stellar mass, 
to minimise the assumptions made about stellar hosts. For Class I 
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Figure 7. As Fig. 6 for detected (top) and non-detected (bottom) Class II 
sources. In addition to the disk source at the central position, nearby disk 
sources, bright enough to appear in the stacked image, are also seen, but are 
clearly separated from the nominal disk position. 




Figure 8 . As Fig. 6 for non-detected Class III sources, with the radial pro¬ 
files shown separately for each wavelength. Bright nearby sources have 
been removed from the stacked data (see text). 
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Figure 9. Disk mass (at 20 K) as a function of spectral type for Class I (open 
squares), Class II (closed squares) and Class III (open circles). Red crosses 
mark known multiple sources. 


Figure 10. Disk mass (at 20 K) calculated from 850 /im flux against that 
calculated from 450 /im flux. Disk class and multiple sources as in Fig. 9, 
with magenta circles and cyan triangles marking disks of spectral type M 
and K, respectively. The dotted line shows the 1:1 ratio. 


sources with a range of possible spectral types, the median value 
has been used. 

Fig. 9 shows the disk mass, calculated assuming a dust tem¬ 
perature of 20 K, as a function of spectral type. Since binary or 
higher multiple disk sources may have multiple and/or circumbi- 
nary disks, we identify known binaries in the plot. The multi¬ 
ple sources all have late spectral types, and with the exception of 
MHO 1/2, are at the lower end of the mass distribution. For both 
Classes, M-type stellar hosts have a wider distribution of disk mass 
than K-type stellar hosts, with stellar hosts of spectral types M2 
to M6 having the largest range of mass. For the Class II sources, 
with the exception of the anomalous CY Tau and multiple MHO 1/2 
sources, M-type stellar hosts have lower masses than K-type stel¬ 
lar hosts. These data support results from studies of late spectral 
type PMS stars in Taurus at high resolution, which have similarly 
found evidence that disk masses decrease over the limited spectral 
range of K-M stars, with a difference in flux observed at the M2- 
M3 spectral type (Schaefer et al. 2009) and Mj oc M* in a linear 
scaling (Andrews et al. 2013). 

Fig. 9 shows a slight trend in the opposite direction for the 
Class I disk sources, with later spectral types extending to larger 
masses, which may be associated with a difference in evolutionary 
timescales for more massive stellar hosts. Given the difficulty in 
determining spectral types for Class I sources, however, this trend 
is much less clearly defined than for the Class II sources. 

Fig. 10 shows the disk mass derived at each wavelength, 
calculated assuming T ( j us t=20K, with source Class and multiple 
sources marked as in Fig. 9. Additionally, magenta circles and cyan 
triangles mark disks with stellar hosts of spectral type M and K, re¬ 
spectively. Disks with stellar hosts of spectral type K consistently 
have higher 450 /im disk masses than 850 /im disk masses. For 
disks with stellar hosts of spectral type M, and for the different 
disk classes, there are no clear trends in the mass ratio. Although 
variations in dust opacity are the largest source of uncertainty in 
disk mass calculations, these results suggest that there could be an 
intrinsic difference in the dust opacity, or the dust temperature, or 
a combination of both, for disks surrounding K-type stellar hosts 
compared to those surrounding M-type stellar hosts. Earlier type 
stars have higher luminosities, which generally leads to higher dust 
temperatures in the disk (Andrews et al. 2013). Once complete, the 


full GBS dataset will provide a much larger sample with which to 
investigate this difference. 


5 CONCLUSIONS 

• CO contamination in the continuum at 850 /im is negligible in 
L1495, with only five sources showing contamination levels over 
5%, all of which are are below 15% 

• Across the L 1495 region, the spectral slope a = 2.8 ±0.6, 
suggesting that much of the dust detected by SCUBA-2 has un¬ 
dergone grain growth. 

• Class I and Class II disks have similar distributions of a val¬ 
ues, with four sources having 2.5 < a < 3, and eight sources with 
1.4 < a < 2.5. These results imply that j8 ~1 for the majority of 
disk sources, but that both Classes show evidence either for small, 
optically thick disks, or rapid grain growth. 

• Class I sources have slightly higher dust colour tempera¬ 
tures than Class II sources, although both Classes have similar 
ranges of dust temperature, 7-40 K. In two Class I sources (IRAS 
04108+2803 A and IRAS 04108+2803 B) the 450/850 gm flux ra¬ 
tio is approaching the upper limit of 50 K, where the dust colour 
temperature is unconstrained, particularly for low j 3 values. 

• In L1495, 2/29 Class II sources have disk masses 
>20Mjup. If masses are calculated from the dust colour tem¬ 
perature, then the two coldest sources have much larger masses i.e., 
3/29 Class II sources and 1/8 Class I sources have Mj >20Mjup. 

• We detect a higher percentage of disks with stellar hosts of 
spectral type K than those of spectral type M, with detection rates 
of 50% and 33% respectively at 850 /im. At 450^1 m, we detect the 
same 50% of objects with stellar hosts of spectral type K, but only 
16% of those of spectral type M. Variations in the dust opacity lead 
to the largest uncertainties in disk mass calculations, but our results 
suggest there is an intrinsic difference in the dust opacity, and/or 
the dust opacity between disks with stellar hosts of spectral type 
M and K. The latter all have higher disk masses calculated using 
450 /im fluxes than calculated using 850 /dm data. 

The Legacy datasets forming the JCMT Gould Belt Survey 
will be used to extend this analysis across a range of star-forming 
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regions, so that comparisons between regions and star-forming en¬ 
vironments can be made. 

The largest uncertainty in calculating disk masses is in de¬ 
termining the dust opacity, particularly for observations where the 
disks are not resolved. Research in this field will make significant 
progress when the full capabilities of ALMA become available, 
and nearby disks can be imaged at millimetre/submillimetre wave¬ 
lengths on scales of a few AU. 
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Figure Al. 850 /im (left) and 450 /im (right) dust continuum towards 
J04153566+2847417 Class I 



Figure A4. as Fig. Al IRAS04111+2800G Class I 



Right ascension (J2000) Right ascension (J2000) 


Figure A2. as Fig. Al SSTTau041831.2+28261 Class I 


Figure A5. as Fig. Al IRAS04108+2803A Class I 


APPENDIX A: DISK CONTINUUM IMAGES 

Figures A1-A56 show thumbnail images of the dust continuum at 
the nominal position for each disk at 850 /im and 450 /im. The flux 
is shown in mJy pixel -1 , and circles show the 20 /r aperture used to 
measure fluxes. The centrally-aligned source is the one named in 
the caption. 





Figure A7. as Fig. Al CoKuTau/1 Class I 



Figure A3, as Fig. Al IRAS04108+2803B Class I 
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Figure A9. as Fig. A1 V410X-ray7 Class II 


Figure A14. as Fig. A1 CWTau Class II 




Figure A10. as Fig. A1 J04141 188+2811535 Class II 


Figure A15. as Fig. A1 CIDA1 Class II 





Right ascension (J2000) 


Figure All. as Fig. A1 V773TauA+B Class II 


Figure A16. as Fig. A1 MHO 1/2 Class II. A larger 24" aperture was used 
for this source. 



Figure A12. as Fig. A1 FMTau Class II 


Figure A17. as Fig. A1 MH03 Class II 






Figure A13. as Fig. A1 FNTau Class II 


Figure A18. as Fig. A1 FOTauA+B Class II 


































































































14 The GBS Taurus Team 



Figure A19. as Fig. A1 J04153916+2818586 Class II 


Figure A24. as Fig. A1 KPNOIO Class II 



Figure A20. as Fig. A1 J04155799+2746175 Class II 


Figure A25. as Fig. A1 V41 OX-ray 1 Class II 



Figure A21. as Fig. A1 J04161210+2756385 Class II 


Figure A26. as Fig. A1 V410Anonl3 Class II 





Figure A22. as Fig. A1 J04163911+2858491 Class II 


Figure A27. as Fig. A1 DDTauA+B Class II 



Figure A23. as Fig. A1 CYTau Class II 


Figure A28. as Fig. A1 CZTauA+B Class II 
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Figure A29. as Fig. A1 V410X-ray2 Class II 


Figure A34. as Fig. A1 FQTauA+B Class II 



44s 43s 42s 41s 40s 39s4hl8m38s 

Right ascension (J2000) 



Figure A30. as Fig. A1 V892Tau Class II 




Figure A31. as Fig. A1 LR1 Class II 


Figure A36. as Fig. A1 FRTau Class II 




Figure A32. as Fig. A1 V410X-ray6 Class II 


Figure A37. as Fig. A1 J04201611+2821325 Class II 



Figure A33. as Fig. A1 KPNQ12 Class II 
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Figure A38. as Fig. A1 LkCal Class III 


Figure A43. as Fig. A1 KPNOl Class III 



Figure A39. as Fig. A1 Anonl Class III 


Figure A44. as Fig. A1 J04161885+2752155 Class III 



Figure A40. as Fig. A1 XEST20-066 Class III 


Figure A45. as Fig. A1 LkCa4 Class III 



55s 54s 53s 52s 51s 50s4hl4m49s 

Right ascension Q2000) 




Figure A41. as Fig. A1 XEST20-071 Class III 


Figure A46. as Fig. A1 LkCa5 Class III 



Figure A42. as Fig. A1 CIDA2 Class III 


Figure A47. as Fig. A1 V410Xray3 Class III 
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Figure A48. as Fig. A1 V410Anon25 Class III 


Figure A53. as Fig. A1 KPNQ2 Class III 



Figure A49. as Fig. A1 KPNOl 1 Class III 


Figure A54. as Fig. A1 V410X-ray5a Class III 




Figure A50. as Fig. A1 V410TauA+B+C Class III 


Figure A55. as Fig. A1 J04144158+2809583 Class III 



44s 43s 42s 41s 40s 39s 38s4hl8m 7s 43s 42s 41s 40s 39s 38s4hl8m37s 

Right ascension (J2000) Right ascension (J2000) 



Figure A51. as Fig. A1 V410X-ray4 Class III 


Figure A56. as Fig. A1 J04161726+2817128 Class III 



Figure A52. as Fig. A1 Hubble4 Class III 








































































































